DISCOVERY OF THE B c MESON 
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We report on the first observation of the bottom-charm mesons B c through 
the decay mode B^r — > J/ipl v in 1.8 TeV pp collisions using the CDF detector 
at the Fermilab Tevatron. A fit of background and signal contributions to the 
observed J/ip£ mass distribution yielded 20.4±5g events from B c mesons. A fit to 
the same distribution with background alone was rejected at the level of 4.8 standard 
deviations. We measured the B+ mass to be 6.40 ± 0.39 ± 0.13 GeV/c 2 and the 
B+ lifetime to be 0.46±[^ (stat.) ± 0.03 (syst.) ps. The measured production cross 
section times branching ratio for B~^~ — > J '/vjj £ + v relative to that for B + — > J /ip K + 
is 0.132±g;g|^ (stat.) ± 0.031 (syst.) ±8:820 (lifetime). 



1 Introduction 



The .B+ meson is the lowest-mass bound state of quarkonium states containing a charm 
quark and a bottom anti-quark |IJ. Since this pseudoscalar ground state has non-zero 
flavor, it has no strong or electromagnetic decay channels. It is the last such meson 
predicted by the Standard Model. It decays weakly yielding a large branching fraction to 
final states containing a J/ip 0, |3|, §], ||. Non-relativistic potential models predict a B c 
mass in the range 6.2-6.3 GeV/c 2 || [7|. In these models, the c and b are tightly bound 
in a very compact system and have a rich spectroscopy of excited states. 

The production of B c mesons has been calculated in perturbative QCD. At transverse 
momenta pt large compared to the B c mass the dominant process is that in which a b is 
produced by gluon fusion in the hard collision and fragmentation provides the c. At lower 
Pt a full calculation || shows that the dominant process is one in which both the 
b and c quarks are produced in the hard scattering. These and other calculations M, 0, 



lOl pi] , |12il provide inclusive production cross sections along with distributions in px and 



other kinematic variables. 

We expect three major contributions to the B c decay width: b —>■ cW + with the c as a 
spectator, leading to final states like J/ipn or J/ipiu] c — > sW + , with the b as spectator, 
leading to final states like B s n or B s £u; and cb —>■ W + annihilation, leading to final 
states like D K , r v T or multiple pions. Since these processes lead to different final states, 
their amplitudes do not interfere. When phase space and other effects are included, the 
predicted lifetime is in the range 0.4-1.4 ps |2|, [13], [14], [15], [16], [H]]. Because of the wide 



range of predictions, a B c lifetime measurement is a test of the different assumptions made 
in the various calculations. Several authors have also calculated the B c partial widths to 
semileptonic final states f|, ||, [|, |^, [18|] . 

Limits on the B c production have been placed by various experimental searches at LEP 



19| , [20| , ^TJ . A prior CDF search placed a limit on B c production in the B^ — > J/ipir' 



decay mode |22 



We report here on the first observation of B c mesons produced in 1.8 TeV pp collisions 
at the Fermilab Tevatron collider using a 110 pb _1 data sample collected with the CDF 
detector. A more detailed description of this work can be found in Ref. |[23|| . We searched 
for the decay channels B+ — > J/ip^f and £>+ — > J/ipe + v followed by J/ip — ■> 
A Monte Carlo calculation of B c production and decay to J/ipiu showed that, for an 
assumed B c mass of 6.27 GeV/c 2 , 93% of the J/ip £ final state particles would have J/ip£ 
masses with 4.0 < M(J/ip£) < 6.0 GeV/c 2 . We refer to this as the signal region, but we 



accepted candidates with M(J/ifj£) between 3.35 and 11 GeV/c 2 . 

We have described the CDF detector in detail elsewhere |24|, p5| . The tracking system 
for CDF gives a transverse momentum resolution Spt/pt = [(0.0009 xp T ) 2 + (0.0066) 2 ] 1 / 2 , 
where pt is in units of GeV/c. The average track impact parameter resolution relative 
to the beam axis is (13 + (40/pr)) fim in the plane transverse to the beam [26]. An 
online di-muon trigger and subsequent offline selection yielded a sample of about 196,000 
J/ip — > mesons. 



2 Event Selection 

We searched for the B c through — > J/ip £ + v decays. These decays have a very simple 
topology: a decay point for J/ip — > [i + fi~ displaced from the primary interaction point 
and a third track emerging from the same decay point. This J/ip + track sample included 
B£ — > J/ipe + i>, B+ — > J/tfjfi + v, B + — > J/ipK + , and background from various sources. 
We subjected the three tracks to a fit that constrained the two muons to the J/ip mass 
and that constrained all three tracks to originate from a common point. A measure of 
the time between production and decay of a B c candidate is the quantity ct*, defined as 

M(J/ip£)-L xy (J/ip£) 

\pt(J/M\ [) 

where L xy is the distance between the beam centroid and the decay point of the B c can- 
didate in the plane perpendicular to the beam direction and projected along the direction 
of the J/ip£ combination in that plane, M(J/ijj£) is the mass of the tri-lepton system, 
and Pt{J/iI> £) is its momentum transverse to the beam. Our average uncertainty in the 
measurement of ct* is 25 /im. We required ct* > 60 /im. 

B + — > J/ip K + candidates were identified by a peak in the fi + [i~K + mass distribution 
centered at M(B + ) = 5.279 GeV/c 2 with an r.m.s. width of 14 MeV/c 2 . (See Fig. 2 of 
Ref. pBj.) The peak contained 290 ± 19 events after correction for background. Events 
within 50 MeV/c 2 of M(B + ) were eliminated as B+ — > J/ip£ + v candidates. 

Electrons were identified by the association of a charged-particle track with p? > 
2GeV/c and an electromagnetic shower in the calorimeter. Additional information for 
identifying electrons was obtained from specific ionization in the tracking chambers and 
from the shower profile in proportional chambers embedded in the electromagnetic calorime- 
ter. Muons from J/ip decay were identified by matching a charged-particle track with 
Pt > 2 GeV/c to a track segment in muon drift chambers outside the central calorimeter 



(5 to 9 interaction lengths thick depending on angle). The third muon was required to 
have a transverse momentum exceeding 3 GeV/c and to pass through an additional three 
interaction lengths of steel to produce a track segment in a second set of drift chambers. 
We found 23 — ► J/ip e + v candidates of which 19 were in the signal region, and we 
found 14 B^ — > J/ip /i + v candidates of which 12 were in the signal region. 



3 Background Determination 

The major contributions to backgrounds in the sample of B c candidates come from 
misidentification of hadron tracks as leptons {i.e. false leptons) and from random com- 
binations of real leptons with J/ip mesons. There are three significant sources of false 
leptons: hadrons that reach the muon detectors without being absorbed; hadrons that 
decay in flight into a muon in advance of entering the muon detectors; and hadrons that 
are falsely identified as electrons. In one type of random combination, electrons from 
photons that convert to e + e~ pairs in the material around the beam line or from Dalitz 
decay of tt° contribute to a "conversion background" when the other member of the pair 
remains undetected. The other type of random combination involves a B that has decayed 
into a J/ip and an associated B that has decayed semileptonically (or through semilep- 
tonic decays of its daughter hadrons) into a muon or an electron. The displaced J/ip 
and the lepton can accidentally appear to originate from a common point. A number 
of other backgrounds |23j were found to be negligible. From a combination of data and 
Monte Carlo calculations, we determined the J/ip + track mass distribution for each of 
the sources of background. As a check of our background calculations, we verified that 
we are able to predict the number of events and mass distribution in an independent, 
background-rich sample of same-charge, low-mass lepton pairs. (See Fig. 27 in Ref. [p3|j.) 

The topology for candidate events was verified by applying all selection criteria except 
the requirement that the third track intersect the J/ip vertex. The impact parameter 
distribution between the third track and the J/ip vertex has a prominent peak at zero, 
demonstrating that, for most candidate events, the three tracks arise from a common 
vertex. (See Fig. 28 in Ref. f23[.) 

In Table [L] we summarize the results of the background calculation and of a simul- 
taneous fit for the muon and electron channels to the mass spectrum over the region 
between 3.35 and 11 GeV/c 2 . Figure [I] shows the mass spectra for the combined J/ipe 
and J /ip /i candidate samples, the combined backgrounds and the fitted contribution from 
B+ — > J /ipi + v decay. The fitted number of B c events is 20 A~^'^. 



Table 1: B c Signal and Background Summary 





3.35 < M(J/iP£) < 11.0 GeV/c 2 




J /ip e Events 


J /ip H Events 


False Electrons 


4.2 ±0.4 




Undetected Conversions 


2.1 ± 1.7 




False Muons 




11.4 ±2.4 


BB bkg. 


2.3 ±0.9 


1.44 ±0.25 


Total Background (predicted) 


8.6 ±2.0 


12.8 ±2.4 


(from fit) 


9.2 ±2.0 


10.6 ±2.3 


Predicted N(B C -> J/ip e v)/N(B c -> J/ipiv) 


0.5^ 


I ± 0.04 


e and /i Signal (derived from fit) 


12.010 


8-4±|;I 


Total Signal (fitted parameter) 


20.4l|i 


Signal + Background 


21.2 ±4.3 


19.0 ±3.5 


Candidates 


23 


14 


Probability for background 






alone to fluctuate to the 


0.63 x 10~ 6 


apparent signal of 20.4 events 







4 Test of Significance and Determination of the B c 
Lifetime and Mass 

To test the significance of this result, we generated a number of Monte Carlo trials with the 
statistical properties of the backgrounds, but with no contribution from B c mesons. These 
were subjected to the same fitting procedure to determine contributions consistent with 
the signal distribution arising from background fluctuations. The probability of obtaining 
a yield of 20.4 events or more is 0.63 x 10 -6 , equivalent to a 4.8 standard-deviation effect. 

To check the stability of the B c signal, we varied the value assumed for the B c mass. 
We generated signal templates, i.e. Monte Carlo samples of B+ — > J/ip £ + v, with various 
values of M(B C ) from 5.52 to 7.52 GeV/c 2 . The signal template for each value of M(B C ) 
and the background mass distributions were used to fit the mass spectrum for the data. 
This study established that the magnitude of the B c signal is stable over the range of 
theoretical predictions for M(B C ), and the dependence of the log-likelihood function on 
mass yielded M(B C ) = 6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . 

We obtained the mean proper decay length cr and hence the lifetime r of the B c meson 
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Figure 1 : Histogram of the J/ip£ mass that compares the signal and background contri- 
butions determined in the likelihood fit to the combined data for J/ip e and J / ip fi. Note 
that the mass bins vary in width. The total B c contribution is 20. 4^5 events. The inset 
shows the behavior of the log-likelihood function —2 ln(L) vs. the number of B c mesons. 



from the distribution of ct*. We used only events with 4.0 < M(J/ip£) < 6.0 GeV/c 2 , 
and we changed the threshold requirement on ct* from ct* > 60 /iin to ct* > —100 /mi for 
this lifetime measurement. This yielded a sample of 71 events, 42 J/ijje and 29 J/ipfi. 
We determined a functional form for the shapes in ct* for each of the backgrounds. To 
these, we added a resolution-smeared exponential decay distribution for a B c contribution, 
parameterized by its mean decay length ct. Because of the missing neutrino, the proper 
decay length ct for each event differs from ct* of Eq. |T[ We convoluted the exponential in ct 
with the distribution of ct* jet derived from Monte Carlo studies. Finally, we incorporated 
the data from each of the candidate events in an unbinned likelihood fit to determine the 
best-fit value of ct. The data and the signal and background distributions are shown in 
Fig. H, and the result is: 

CT = 137±49 ( stat -) ± 9 ( s Y st ) V m ( 2 ) 



0.461°;^ (stat.) ± 0.03 (syst.) 



ps 



(3) 




Figure 2: The distribution in ct* for the combined J/ipfi and J/ipe data along with the 
fitted curve and contributions to it from signal and background. The inset shows the 
log-likelihood function vs. ct for the B c . 



From the 20.4 B c events and the 290 B + — » J/ip K + events, we calculated the B c pro- 
duction cross section times the B+ — > J/ip £ + v branching fraction, a-BR(B^ — > J/ip 
relative to that for the topologically similar decay B + — > J/ip K + . Systematic uncertain- 
ties arising from the luminosity, from the J / ip trigger efficiency, and from the track-finding 
efficiencies cancel in the ratio. Our Monte Carlo calculations yielded the values for the 
efficiencies that do not cancel. The detection efficiency for B+ — > J/ip £ + v depends on 
ct because of the requirement that ct* > 60 /xm, and we quote a separate systematic 
uncertainty because of the lifetime uncertainty. We assumed that the branching fraction 
is the same for B+ — > J/ipe + v and B+ — > J/ip fi + v. We multiply the 20.4 events by a 
factor 0.85 ± 0.15 to correct for other B c decay channels such as B c — ► ip' £v [p3[]. We find 



n ^ e ^ s T{b) : Tr^Zj'/Ik) = ° 132 ^ (stat) ±0 - 031 (syst - } (lifetime) ' 

(4) 

for B^ and B + with transverse momenta px > 6.0 GeV/c and rapidities \y\ < 1.0. This 
result is consistent with limits from previous searches |L9|, ^(| Figure |3] compares 
phenomenological predictions with our measurements of ct and 1Z( J/ip Iv). Within ex- 
perimental and theoretical uncertainties, they are consistent. 
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Figure 3: The point with 1-standard-deviation contour shows our measured value of the 
a ■ BR ratio plotted at the value we measure for the B c lifetime. The shaded region 
represents theoretical predictions and their uncertainty corridors for two different values 
of the semileptonic width T s j. = T(B C — > JjijjEu) based on Refs. and j|. The other 
numbers assumed in the theoretical predictions are V c b = 0.041 ±0. 005 |?7), a(B+) / a(b) = 
1.3 x 10- 3 [@], ^ff- = 0.378 ±0.022 |27|, BR{B+ -> J/^K+) = (1.01 ±0.14) x lO" 3 0. 



5 Conclusions 

In conclusion, we report the observation of B c mesons through their semileptonic decay 
modes, B c — ► J/ip £X where I is either an electron or a muon. We measured the B c mass 
and the product of its production cross section times semileptonic branching fraction, 
which confirm phenomenological expectations. We measured a B c lifetime consistent with 
calculations in which the decay width is dominated by the decay of the charm quark. 
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